Technique for Direct eV-Scale Measurements of the Mu and Tau Neutrino Masses 

Using Supernova Neutrinos 
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Early black hole formation in a core-collapse supernova will abruptly truncate the neutrino fluxes. 
The sharp cutoff can be used to make model-independent time-of-flight neutrino mass tests. Assum- 
ing a neutrino luminosity of f 52 erg/ s per flavor at cutoff and a distance of f kpc, SuperKamiokande 
can detect an electron neutrino mass as small as f .8 eV, and the proposed OMN1S detector can 
detect mu and tau neutrino masses as small as 6 eV. This Letter presents the first technique with 
direct sensitivity to eV-scale mu and tau neutrino masses. 
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Introduction: Despite decades of experimental ef- 
fort, the values of the neutrino masses remain elusive. 
While the laboratory bound on the electron neutrino 
mass is about 3 eV [Q , the laboratory bounds on the mu 
and tau neutrino masses are much weaker: 170 keV ||] 
and 18 MeV H , respectively. Only recently have neutrino 
oscillation experiments found strong evidence for nonzero 
differences of squared neutrino masses. Once discovered, 
the values of the neutrino masses may provide important 
clues to physics beyond the Standard Model. In some sce- 
narios, e.g., with the see-saw mechanism j|, the mu and 
tau neutrino masses are expected to be much larger than 
the electron neutrino mass. If they are at the eV scale or 
greater, the neutrino masses could also be important cos- 
mologically as a component of the long-sought dark mat- 
ter. It is therefore crucial to devise direct tests of the mu 
and tau neutrino masses with sensitivity reaching the eV 
scale. While neutrino mass tests based on cosmological 
considerations may reach the eV scale, they are indirect 
(no neutrinos are detected) and depend upon the other 
cosmological parameters being independently known || . 

The best known possibility for directly measuring the 
mu and tau neutrino masses is by time-of-flight measure- 
ments of supernova neutrinos, comparing the arrival time 
of the mu and tau neutrinos to that of the electron neu- 
trinos. However, this is complicated by the long intrinsic 
duration (~ 10 s) of the neutrino signal and the fact that 
its detailed characteristics are mo del- dependent. Bea- 
com and Vogel have shown that a technique based on 
the average arrival times (t) is model-independent and is 
sensitive to delays as small as ~ 0.1 s ||. This would 
allow detection of mu or tau neutrino masses down to 
45 eV in SuperKamiokande (SK) and 30 eV in the Sud- 
bury Neutrino Observatory (SNO). If the mu and tau 
neutrino masses (strictly speaking, those of the relevant 



mass eigenstates) are nearly degenerate, as suggested by 
the atmospheric neutrino results 0, then the sensitiv- 
ity would improve by about a/2- Unfortunately, it seems 
difficult to improve the results with this technique, since 
the mass sensitivity scales with the detector mass Md 
as l/Mp 4 ||. To reach the few-eV scale would require 
detectors 10 4 times larger, which seems impossible. 

In this Letter, we discuss a new time-of-flight technique 
for measuring neutrino masses that can reach the eV 
scale. This technique is applicable if the proto-neutron 
star forms a black hole early enough to abruptly termi- 
nate the neutrino signal. We state only our most impor- 
tant results; the details will be discussed at length in a 
forthcoming paper ||. 

Expected Neutrino Signal: We consider black 
hole formation which occurs soon (~ 1 s) after core col- 
lapse (other scenarios are considered in Ref. ||). Black 
hole formation is triggered by accretion, which drives the 
proto-neutron star mass above the maximum stable neu- 
tron star mass. The neutrino signal expected in this sce- 
nario has been studied by Burrows Q and Mezzacappa 
and Bruenn |1C( |. In these models, the neutrino lumi- 
nosities were fairly constant at more than 10 52 erg/s per 
flavor until abruptly terminated by black hole formation. 
In fact, the transition should have a nonzero duration, of 
order the light crossing time 2R/c ~ 0.1 ms, as the proto- 
neutron star radius shrinks to that of the final black hole. 
During the transition, the gravitational redshift, origi- 
nally ~ 10%, rapidly diverges, truncating the neutrino 
signal. Using a singularity-avoiding code, Baumgarte et 
al. [jllfl studied the transition and found its duration to 
be 0.5 ms. Thus, we can consider the neutrino fluxes to 
be sharply and simultaneously terminated. 

The results below assume a luminosity Lbh = 10 52 
erg/s per flavor at the cutoff time tsH, and a distance 
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D — 10 kpc. We assume the following temperatures: 
T = 3.5 MeV for u e , T = 5 MeV for v e , and T = 8 MeV 
for v^, v T , v^, and v r . It will be shown that the necessary 
quantities can be measured in a realistic situation. 

Neutrino Mass Effects: At lowest order, a neutrino 
with mass m (in eV) and energy E (in MeV) will have 
an energy-dependent delay (in s) relative to a massless 
neutrino in traveling over a distance D (in 10 kpc): 



(1) 



At(E) = 0.515 (— ] D. 



The distance is scaled by the approximate distance to 
the Galactic center, though a supernova may be detected 
from anywhere in the Galaxy. For the smallest detectable 
masses, the delay effects will be visible only after the 
sharp cutoff, where no events are otherwise expected. 
Since the delays are very small, the luminosities and tem- 
peratures can be taken as constant over the short interval 
before t BH . The event rate for t > t BH is M: 



dN 
~dt 



(t) = C 



L 



BH 



10 51 erg/s 



dEf(E) 



a(E) 



10" 



(2) 



where f(E) is the neutrino energy spectrum and cr{E) 
the cross section. The upper limit E max on the integral 
allows only delays as large as t — tsn, i.e., 
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where the units are as in Eq. (Q). Note that the time and 
neutrino mass dependence appear only through E max . 



For t < t B H, E„ 



oo, and the rate is constant. If the 



neutrino energy can be measured, as for some charged- 
current reactions, then the event rates for different neu- 
trino energies can easily be obtained. For an H2O detec- 
tor, the constant C is 



C H 20 = (1.74/s) 



Mr 



1 kton 



10 kpc 
D 



2 r 



1 MeV 



(E) 



(4) 



For a Fermi-Dirac spectrum, (E) = 3.15T. The constant 
for a 208 Pb detector can be obtained by scaling by the 
relative number of targets/kton, i.e., 18/208. 

The expected number of delayed counts after t B H can 
be calculated using Eq. (ph. This will be useful when tsn 
can be measured independently. It can be shown § that 
this has the very simple form: 



dN ( m 

N dei = -j-(t BH ) x 0.515 — 

dt \ En 



D. 



(5) 



where the event rate is in s _1 , and the other units are 
as in Eq. (Q). This formula would obviously be true if 
only a single energy contributed and the sharp cutoff in 
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FIG. 1. The event rates due to v e + p — > e + + n in 
SK, for different ranges of the neutrino energy: "Low" 
(0 < E < 11.3 MeV, contains 2.4 events past the true 
tan), "Mid" (11.3 < E < 30 MeV, 4.8 events), "High" 
(30 < E < 00 MeV, 0.5 events), and "All" (all energies, 7.7 
events). Note that only the rate after about tsn is shown, 
and that the range of t — tsH is very short. 



the event rate were rigidly translated by the delay. But 
it is remarkable and very convenient that it is still true 
even when there is a spectrum of energies and the event 
rate develops a decaying tail past the cutoff (as in Figs. |l| 
and|^). The physical significance of the "central" energy 
E c is that it is (to an excellent approximation) simply 
the Gamow peak of the falling thermal spectrum and the 
rising cross section. As derived, this is an exact result. 

Electron Neutrino Mass: We first consider the 
measurement of tBH and m Uc using the v e + p — > e + +n 
events in the 32-kton SK detector. For T = 5 MeV, the 
thermally-averaged cross section (for the sum of the two 
protons in H 2 0) is 44. x 10~ 42 cm 2 [Q. The event rate 
at or before tsn is thus ~ 1500 s _1 . After tBH, the rate 
is zero if m„ c = and will develop a tail if ro„ t > 0. 

For a sharp edge, the edge position can be determined 
with an error given by the reciprocal of the event rate 
before the edge, i.e., the event spacing ||,|l3|]. If we knew 
that m Ue — 0, then tBH would be determined to ~ 1 
ms. More realistically, a mass as large as the laboratory 
bound, m Ve < 3 eV Q, would cause delays as large as 40 
ms, so that the extracted tsn would be too large. 

However, we can simultaneously measure m„ e and tBH 
by splitting the v e +P — * e + +n data into different ranges 
of neutrino energy (using E v ~ E e + 1.3 MeV). These are 
defined in the caption of Fig. |l|. The Low group must be 
excluded from consideration because these events have 
positron total energy less than 10 MeV, and can be con- 
fused with the 5-10 MeV gammas from neutral-current 
reactions on 16 O |1^] . The High group has very little 
delay and will thus primarily be sensitive to tsn- Then 
the Mid group will determine m„ e , by counting events 
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FIG. 2. The results for the combined 1-n neutral-current 
event rate due to , v T , v^, and v r in OMNIS. Note that only 
the rate after about tsH is shown. Before tBH, other reac- 
tions contribute about 20% of the total neutron rate; they are 
not included here, and will have to be statistically subtracted 
from the measured rate. The mu and tau neutrino masses are 
assumed degenerate . The m = case is drawn with a solid 
line. The m = 6.1 eV case, with 2.3 events expected in the 
tail, is the first case that can be reliably distinguishable from 
m = 0, and is drawn with a long-dashed line. The results for 
other masses are drawn with dotted lines. 



delayed past the tsH determined by the High group. 

In Fig. |l|, we show a possible analysis for the case of 
to„ c =1.8 eV. In the High group, the number of events in 
the tail is < 1, so the cutoff appears sharp and is specified 
to within ~ 2 ms. This uncertainty affects the expected 
number in the Mid group by ~ 2 events. Even so, one can 
still reliably see a few delayed counts after the measured 
tBH , enough to establish a nonzero mass (the statistics 
are discussed in more detail below). A more sophisticated 
fit would improve our results somewhat, and we assume 
a final uncertainty on tsn of about I ms. For a super- 
nova in which the neutrino fluxes are not truncated by 
black hole formation, SK could detect an electron neu- 
trino mass as small as ~ 3 eV |ll| . 

Mu and Tau Neutrino Masses: We consider mu 
and tau neutrino detection in OMNIS, a proposed super- 
nova neutrino detector based on lead and iron [jl6| . Since 
their energies are below the charged-current thresholds, 
supernova mu and tau neutrinos can be detected only by 
their neutral-current interactions. On the other hand, 
due to the temperature hierarchy, they will dominate the 
neutral-current yields. In OMNIS, the dominant neutral- 
current reaction is the spallation of single neutrons from 
lead. The neutrons could be detected by capture in a 
gadolinium-doped liquid scintillator, which yields an 8- 
MeV gamma cascade in about 0.030 ms (much smaller 
than typical mass delays). 

For T — 8 MeV, the thermally-averaged cross section 
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FIG. 3. The expected number of delayed counts N,iei in 
OMNIS as a function of the neutrino mass. The points are 
obtained by direct numerical integration of Eq. (Q) . The "+" 
indicates the smallest discernible mass at the 90% CL. The 
solid line is obtained with Eq. (^|), using E c = 40.7 MeV, the 
Gamow peak energy. 



for the sum of and (or v T and v T ) on 208 Pb, includ- 
ing the I-neutron spallation probability, is 760 x 10~ 42 
cm 2 [|l7j. The cross sections on 206 Pb and 207 Pb, which 
together comprise 46% of natural lead, are expected to 
be similar For a supernova at 10 kpc in which the 
neutrino fluxes are not cut off by black hole formation, we 
assume that OMNIS will have ~ 1000 1-neutron neutral- 



current events due to v u 



and v T on lead (the 



events on iron are not included in our calculations) . This 
goal could be met with a 2.2 kton lead detector with per- 
fect neutron detection efficiency. A realistic design based 
on 4 kton of lead and 10 kton of iron, and with about 
this many events, is described by Boyd fic| . 

In Fig. ||, we plot the relevant neutral-current rate for 
different values of the neutrino mass, calculated using 
Eq. (^). In Fig. ||, we plot the number of delayed events 
Ndei as a function of the neutrino mass, using Eq. (|^) 
and by direct integration. Equation (|J) is remarkable for 
its simplicity, and also because it is written in terms of 
measurable quantities. The cutoff time tBH will be mea- 
sured in SK. The neutral-current event rate at or before 
tsH will be measured in OMNIS, as will Ndei- The cen- 
tral energy E c depends on the mu and tau neutrino tem- 
perature, which can be estimated by the neutral-current 
yields on different targets ||. We assume that the dis- 
tance D can be determined by astronomical means. 

Given the measured value of Ndei , Eq. (|^) can be imme- 
diately solved for the best-fit neutrino mass. If Ndei = 
is measured, then the best-fit mass is m = 0, and an 
upper limit can be placed. An expectation of 2.3 counts 
fluctuates down to counts only 10% of the time. Thus, 
setting Ndei — 2.3, an upper limit on the mass mu m is 
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obtained. This is the largest mass, given the expected 
Poisson statistics, that could be confused with the mass- 
less case. For the present case, this is 6.1 eV. 

Since the fractional error on Ndei due to Poisson statis- 
tics is large (~ l/\^3 ~ 65%), errors on other inputs 
are expected to be irrelevant. The uncertainty on tsn 
from SK is assumed to be about 1 ms. From Fig. (g), this 
uncertainty can be seen to change the expected number 
Ndel by ~ 0.2 events, which is negligible. Other possible 
errors, e.g., the detector background, the disregarded 0.5 
ms tail of the luminosity, and v e and v e events after tsH, 
are even less important [pi. 

For a supernova that does not have the sharp cutoff in 
the rate characteristic of black hole formation, the model- 
independent (t) analysis yields an mn m that is inde- 
pendent of the distance D and scales as l/M^ 0. For 
the present case, mu m scales as: 



rnj, 



(E)D 



VeffLBHMD 



(6) 



where <r e ff is the thermally-averaged cross section. In 
terms of absolute sensitivity, these techniques compare 
as 21 eV and 6 eV, respectively. These differences are 
consequences of the sharp cutoff in the neutrino flux. 

Conclusions: If a black hole forms early in a core- 
collapse supernova, then the fluxes of the various flavors 
of neutrinos will be abruptly and simultaneously termi- 
nated when the neutrinospheres are enveloped by the 
event horizon. For a massive neutrino, the cutoff in the 
arrival time will be delayed by At ~ (m/E) 2 relative to 
a massless neutrino. 

The Galactic core-collapse supernova rate is about 
3/century or higher S, and the work of Brown and 
Bethe JL8| suggests that black holes are formed about 
half of the time. In the work of Burrows || and Mez- 
zacappa and Bruenn [[Tot , the neutrino luminosities just 
before black hole formation are very high. These results 
indicate that there is a reasonably good chance that such 
an event could be observed by the present and proposed 
supernova neutrino detectors p]. If so, there are impor- 
tant practical consequences. 

First, since SK can measure the neutrino energy of 
the u e + p — > e + + n events, both tsu and m Ve can be 
measured by the arrival times for different neutrino en- 
ergies. An electron neutrino mass as small as 1.8 eV can 
be detected. Second, although the mu and tau neutrino 
energies are not measured in their neutral-current detec- 
tion reactions, their masses can be measured by counting 
the number of events after tsn ■ In the proposed OMNIS 
detector, a mu and tau neutrino mass (assumed degen- 
erate Q) as small as 6 eV can be detected. This is the 
only known direct technique with eV-scale sensitivity for 
these masses. Third, these results scale with the distance, 
luminosity, and detector mass as yj ' D / ' LbhMd- This fa- 
vorable scaling with the detector mass suggests that it 



would be realistic to consider even larger detectors, in 
order to reach 1 or 2 eV for all three neutrino masses. 
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